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Abstract 

We present an updated measurements of time-dependent CP asymmetries in fully reconstructed 
neutral B decays to several CP eigenstates containing a charmonium meson. The measurements use 
a data sample of (347.5 ± 3.8) x 10 6 T(AS) -> BB decays collected with the BABAR detector at the 
PEP-II B factory between 1999 and 2006. We determine sin2/3 = 0.710 ± 0.034(stat) ± 0.019(syst) 
and |A| = 0.932 ± 0.026(stat) ± 0.017(syst). Both of these results are preliminary. 
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1 INTRODUCTION 



Charge conjugation-parity (CP) violation in the B meson system has been established by the 
BABAR and Belle [2J collaborations. The Standard Model (SM) of electroweak interactions 
describes CP violation as a consequence of a complex phase in the three-generation Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix [3]. In this framework, measurements of CP 
asymmetries in the proper-time distribution of neutral B decays to CP eigenstates containing a 
charmonium and K° meson provide a direct measurement of sin2/3 [3]. The angle f3 is defined as 
arg[— (^.(j^cb)/(^td^tb)]' wnere the Vij are CKM matrix elements. 

In this paper, we report on an updated measurement of sin2/3 in (347.5 ±3.8) x 10 6 T(45) — > BB 
decays collected with the BABAR detector using B° decays to the final states J/ipK®, J/ipK®, 
ip(2S)K®, XciK®, fJ c Kg, and J/ip K*° (K*° — » K^t: ). Since our previous measurement 0, we have 
added a sample of 120 x 10 6 T(4S) — » BB decays, and applied an improved event reconstruction to 
the complete dataset. A new r] c K® event selection has been developed based on the Dalitz structure 
of the rj c — > K®K + tt~ decay. We have also performed a more detailed study of the CP properties 
of our background events resulting in a reduced systematic error. 



2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric 
e + e~ storage ring from 1999 to 2006. This represents a total integrated luminosity of (316.2 ± 
3.5) fb _1 taken at the T(4S) resonance (onpeak), corresponding to a sample of (347.5 ± 3.8) x 
10 6 T(4S) -> BB decays. 

The BABAR detector is described in detail elsewhere @. Charged particles are selected and 
their momenta are measured by a combination of a vertex tracker consisting of five layers of 
double-sided silicon microstrip detectors and a 40-layer central drift chamber, both operating in 
the 1.5 T magnetic field of a superconducting solenoid. We identify photons and electrons using a 
CsI(Tl) electromagnetic calorimeter (EMC). Charged particle identification (PID) is provided by 
an internally reflected ring imaging Cherenkov detector (DIRC) covering the central region of the 
detector, the average energy loss (dE/dx) in the tracking devices, and by the EMC. In addition, 
the instrumented flux return (IFR) containing resistive plate chambers are used for muon and long- 
lived neutral hadron identifications. We use the GEANT4 [7j software to simulate interactions of 
particles traversing the BABAR detector. 



3 ANALYSIS METHOD 



We use information from the other B meson, B tag , in the event to tag the initial flavor of the 
fully reconstructed B candidate. The decay rate /+(/_) for a neutral B meson decaying to a CP 
eigenstate accompanied by a B°(B°) tag can be expressed in terms of a complex parameter A [Sj 



as 



/±(At) 



-|At|/r B o 
4t b o 



(lq=Au>)±(l-2w) x 



21m X 
1 + IAI 2 



sin (Am,j At) 



1+A 



IAI 2 

— L cos (Am d At) 



(1) 



where At = t rec — 
and tagged (B tag/ 



t ta g is the difference between the proper decay times of the reconstructed (B rec ) 
B mesons. t b o is the B° lifetime and Am^ is the mass difference determined 
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from B°-B° oscillations The average mistag probability u describes the effect of incorrect tags, 
and Ao; is the difference between the mistag rate for B° and B . Here we assume that the decay 
width difference AI^ between the neutral B mass eigenstates is zero. The sine term in Eq. ^ is 
due to the interference between the direct decay and the decay after B°-B° oscillation. A non-zero 
cosine term arises from the interference between decay amplitudes with different weak and strong 
phases (direct CP violation) or from CP violation in B°-B° mixing. 

In the SM, CP violation in mixing is negligible, as is direct CP violation for b — ► ccs decays 
that contain a charmonium meson |Hj. Under these assumptions, A = r\fe~ 2% ^ where r/f = ±1 is 
the CP eigenvalue of the final state /. Thus, the time-dependent CP-violating asymmetry is 

Acp{At) = u0) + /"(At) a ~ (1 ~ 2u)r]f s ' m2(3s ' m {Amd At) - (2) 

We reconstruct a sample of neutral B mesons, Bqp, decaying to the rjf = — 1 final states ^U] 
J/tpK% ip(2S)K°, XciK° 3 and r] c K% and the rjf = +1 final state J^K®. The J/tp and ip(2S) mesons 
are reconstructed through their decays to e + e~ and /x + fi~; the ip{2S) is also reconstructed through 
its decay to J/tpir + Tr~ . The Xci meson is reconstructed in the decay mode J/07- The rj c meson 
is reconstructed in the decay mode KgK + 7r~. We also reconstruct the J/i/jK*°(K*° — » K^-k ) 
final state which can be CP-even or CP-odd due to the presence of even (L=0, 2) and odd (£=1) 
orbital angular momenta contributions. Ignoring the angular information in J/ipK*° results in a 
reduction of the measured CP asymmetry by a factor 1 1 — 2R± | , where R± is the fraction of the 
L=l contribution. We have measured R± = 0.233 ± O.OlO(stat) ± 0.005(syst) which gives an 
effective r/f = 0.504 ± 0.033, after acceptance corrections. 

In addition to the CP modes described above, a large sample Pfl av of B° decays to the flavor 
eigenstates D^~h + (h + = ir + ,p + , and af) and J/ip K*°(K*° — > K + n~) is used for calibrating 
the flavor tagging performance and At resolution. We perform studies to measure apparent CP 
violation from unphysical sources using a control sample of B + mesons decaying to the final states 
J/^K ( * )+ , ^(2S)K + , XdK + , and r) c K+. Since the previous BABAR analyses, we apply an improved 
event reconstruction to all the events and we find that on a subset of the data the new and 
old event reconstruction algorithms give consistent results. The event selection and candidate 
reconstruction are unchanged from those described in Refs. with the exceptions described 

below. As in Ref. jS] we reconstruct the B° — > r] c K® and B^ — > rjcK^ modes using only the 
Tf c — > KgK + -K~ decay with 2.91 < rn K o K + 7T - < 3.05GeV/c 2 . We now exploit the fact that the 
r] c decays predominantly through a Kir resonance at around 1430 MeV/ c 2 , and apply the selection 
criteria 1.26GeV/c 2 < m(K%ir-) < 1.63GeV/c 2 and 1.26GeV/c 2 < m(K + ir-) < 1.63GeV/c 2 . 

We calculate the time interval At between the two B decays from the measured separation Az 
between the decay vertices of P rec and B tSug along the collision (z) axis [T^]. The z position of the 
B rec vertex is determined from the charged daughter tracks. The Ptag decay vertex is determined 
by fitting tracks not belonging to the B TCC candidate to a common vertex, employing constraints 
from the beam spot location and the B Tec momentum |12j . Events are accepted if the calculated 
At uncertainty is less than 2.5 ps and |At| is less than 20 ps. The fraction of events satisfying these 
requirements is 95%. The r.m.s. At resolution is 1.1 ps for the 99.7% of events that are not outliers. 

Multivariate algorithms are used to identify signatures of B decays that determine ( "tag" ) the 
flavor of the i?tag at decay to be either a B° or B° candidate. Primary leptons from semileptonic 
B decays are selected from identified electrons and muons as well as isolated energetic tracks. The 
charges of identified kaon candidates are used in a kaon tag. Low momentum pions from D* + 
decays are selected on the basis of their momentum and direction with respect to the thrust axis 
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of -Btag- These algorithms are combined to account for correlations among different sources of 
flavor information and to provide an estimate of the mistag probability for each event. Each event 
whose estimated mistag probability is less than 45% is assigned to one of six tagging categories. 
The Lepton category contains events with an identified lepton; the remaining events are divided 
into the Kaon I, Kaon II, Kaon-Pion, Pion, or Other categories based on the estimated mistag 
probability. For each category i, the tagging efficiency E\ and fraction Wi of events having the 
wrong tag assignment are measured from data (Table ^) . The figure of merit for tagging is the 
effective tagging efficiency Q = J2i ~~ 2roj) 2 = (30.4±0.3) %, where the error shown is statistical 
only. 



Table 1: Efficiencies £j, average mistag fractions Wi, mistag fraction differences Awi = Wi(B°) — 
Wi(B°), and effective tagging efficiency Qi extracted for each tagging category i from the -Bfj av 
sample. 



Category 


£ ' 


(%) 


w 


(%) 


Aw 


(%) 


Q(%) 


Lepton 


8.67 


± 


0.08 


3.0 


±0.3 


-0.2 


±0.6 


7.67 ±0.13 


Kaon I 


10.96 


± 


0.09 


5.3 


±0.4 


-0.6 


±0.7 


8.74 ±0.16 


Kaon II 


17.21 


± 


0.11 


15.5 


±0.4 


-0.4 


±0.7 


8.21 ±0.19 


Kaon-Pion 


13.77 


± 


0.10 


23.5 


±0.5 


-2.4 


±0.8 


3.87 ±0.14 


Pion 


14.38 


± 


0.10 


33.0 


±0.5 


5.2 


±0.8 


1.67 ±0.10 


Other 


9.61 


± 


0.08 


41.9 


±0.6 


4.6 


±0.9 


0.25 ±0.04 


All 


74.60 


± 


0.12 










30.4 ±0.3 



With the exception of the J/ipK® mode, we use the beam-energy substituted mass toes = 
\J (-^beam) 2 ~~ (p*b) 2 *° determine the composition of our final sample, where E^ eam and p* B are the 
beam energy and B momentum in the e + e~ center-of-mass frame. For the J/ipK® mode we use the 
difference AE between the candidate center-of-mass energy and E£ eam . The composition of our 
final sample is shown in Fig. ^ We use events with toes > 5.2GeV/c 2 (AE < 80MeV for J/ifiK®) 
in order to determine the properties of the background contributions. We define a signal region 
5.27 < toes < 5.29GeV/c 2 (\AE\ < lOMeV for J/^K° L ) that contains CP candidate events that 
satisfy the tagging and vertexing requirements as listed in Table [2j 

For all modes except r] c K® and J/ipK® we use simulated events to estimate the fractions of 
events that peak in the toes signal region due to cross-feed from other decay modes (Peaking 
background) . For the r] c Kg mode the cross- feed fraction is determined from a fit to the tukk-k an d 
toes distributions in data. For the J/ip K® decay mode, the sample composition, effective rjf, and 
AE distribution of the individual background sources are determined either from simulation (for 
B — > J/ijjX) or from the m^+£- sidebands in data (for fake J/ip — > l + £~). 

We determine sin2/3 with a simultaneous maximum likelihood fit to the At distribution of the 
tagged Bcp and -Bfl av samples. The At distributions of the Bcp sample are modeled by Eq.^with 
|A| = 1. Those of the -Bfl av sample evolve according to the known frequency for flavor oscillation in 
B° mesons. We assume that the observed amplitudes for the CP asymmetry in the Bqp sample and 
for flavor oscillation in the i?fl av sample are reduced by the same factor 1 — 2w due to flavor mistags. 
The At distributions for the signal are convolved with a resolution function common to both the 
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Figure 1: Distributions for Bcp and i?fl av candidates satisfying the tagging and vertexing require- 
ments: a) toes f° r the final states J/ipK®, ip(2S)Kg, XciKg, and r] c K®, b) AE for the final state 
J/ipK®, c) to es for J/ipK*°(K*° — > KgTT ), and d) toes for the B Rllv sample. In each plot, the 
shaded region is the estimated background contribution. 



Sfi av and Bcp samples, modeled by the sum of three Gaussians J2j. Backgrounds are incorporated 
with an empirical description of their At spectra, containing zero and non-zero lifetime components 
convolved with a resolution function distinct from that of the signal. 

There are 65 free parameters in the fit: sin2/3 (1), the average mistag fractions w and the 
differences Aw between B° and B° mistag fractions for each tagging category (12), parameters 
for the signal At resolution (7), parameters for CP background time dependence (8), and the 
difference between B° and B° reconstruction and tagging efficiencies (7); for -Bfj av background, 
time dependence (3), At resolution (3), and mistag fractions (24). For the CP modes (except for 
J/ipK®), the apparent CP asymmetry of the non-peaking background in each tagging category is 
allowed to be a free parameter in the fit. 

We fix t b o = 1.530 ps, Ato^ = 0.507 ps" 1 9 , |A| = 1, and Ar^ = 0. The determination of the 
mistag fractions and At resolution function parameters for the signal is dominated by the large 
-Bflav sample. We determine background parameters mainly from events outside the peaks in the 
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Table 2: Number of events iVt ag in the signal region after tagging and vertexing requirements, 
signal purity P including the contribution from peaking background, and results of fitting for CP 
asymmetries in the P>cp sample and various subsamples. In addition, results on the -Bfl av and 
charged B control samples test that no artificial CP asymmetry is found where we expect no CP 
violation (sin2/3 = 0). Errors are statistical only. The signal region is 5.27 < th-es < 5.29GeV/c 2 
(|A£| < lOMeV for J/i/)K%). 



Sample 




^tag 


P(%) 


sin2/? 


Full CP sample 




11496 


76 


0.710 ± 0.034 


J/i(>K , iM2S)K»,XciK t i,7 lc 




6028 


92 


0.713 ±0.038 






4323 


55 


0.716 ± 0.080 


J/i>K*°(K*° -> K° s tt°) 




965 


68 


0.526 ± 0.284 


1999-2002 data 




3084 


79 


0.755 ± 0.067 


2003-2004 data 




4850 


77 


0.724 ± 0.052 


2005-2006 data 




3562 


74 


0.663 ± 0.062 


J/^Kl 4>(2S)Kl XdKl 


rj c K% only (77/ = -1) 








J/i/)K° (K s ^ir°ir ) 




4076 


96 


0.715 ± 0.044 




988 


88 


0.581 ±0.105 


i>(2S)K° s vr+vr-) 




622 


83 


0.892 ±0.120 


XdK° s 




279 


89 


0.709 ±0.174 


VcK° s 




243 


75 


0.717 ±0.229 


Lepton category 




703 


97 


0.754 ± 0.068 


Kaon I category 




900 


93 


0.713 ±0.066 


Kaon II category 




1437 


91 


0.711 ±0.075 


Kaon-Pion category 




1107 


89 


0.635 ±0.117 


Pi on category 




1238 


91 


0.587 ±0.175 


Other category 




823 


89 


0.454 ± 0.469 


-Bflav sample 




112878 


83 


0.016 ±0.011 


B + sample 




27775 


93 


0.008 ± 0.017 



rriES an d A.E distribtions, as shown in Fig. ^ 
The fit to the Bqp and B^ v samples yields 

sin2/3 = 0.710 ± 0.034(stat) ± 0.019(syst). 

Figure [2 shows the At distributions and asymmetries in yields between events with B° tags and 
B° tags for the rjf = —1 and rjf = ±1 samples as a function of At, overlaid with the projection 
of the likelihood fit result. We perform a separate fit with only the cleanest rjf = —1 sample, in 
which we treat both |A| and sin2/3 as free parameters. We do not use the modes J/ipK*° and 
J/ipK® to minimize the dependence of the results on the background parametrization. We obtain 
|A| = 0.932 ± 0.026(stat) ± 0.017(syst). The correlation between the coefficients multiplying the 
sin(AmrfAt) and cos(Am^At) terms in Eq. ^is —1.2%. 
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Figure 2: a) Number of rjf = —1 candidates (J/t/jK®, ip(2S)Kg, x&K®, and rj c K®) in the 
signal region with a B° tag (N B o) and with a B° tag (N^ ), and b) the raw asymmetry 
(-/V B o — Ng )/(N B o + iVgo), as functions of At. Figures c) and d) are the corresponding dis- 
tributions for the r\f = +1 mode J/tpK®. All distributions exclude Other-tagged events. The solid 
(dashed) curves represent the fit projections in At for B° (B°) tags. The shaded regions represent 
the estimated background contributions. 



4 SYSTEMATIC UNCERTAINTIES 

The systematic uncertainties on sin2/3 and |A| are summarized in Table El These include the uncer- 
tainties in the level and CP asymmetry of the peaking background, the assumed parameterization 
of the At resolution function, possible differences between the -Bfl av and Bcp tagging performances 
and At resolution functions, knowledge of the event-by-event beam spot position, and the possible 
interference between the suppressed b — > ucd amplitude with the favored b — > cud amplitude for 



13 



some tag-side B decays 14 . In addition, we include the variation due to the assumed values of 
Awid and tb 9 . We also assign the change in the measured sin2/3 as the corresponding systematic 
uncertainties when we let |A| to be a free parameter in the fit and when we set Ar^/r^ = ±0.02, 
the latter being considerably larger than SM estimates JSj. The total systematic error on sin2/3 
(|A|) is 0.019 (0.017). 



Table 3: Systematic uncertainties on sin2/3 and |A|. 



Source 


cj(sin2/3) 




CP backgrounds 


0.007 


0.002 


At resolution function 


0.008 


0.002 


J/ipK® backgrounds 


0.007 


N/A 


Mistag fraction differences 


0.009 


0.007 


Beam spot 


0.008 


0.004 


Am d , t b , Aiyr d , |A| 


0.003 


0.001 


Tag-side interference 


0.002 


0.014 


MC statistics 


0.003 


0.005 


Total systematic error 


0.019 


0.017 



The large Bcp sample allows a number of consistency checks, including separation of the data 
by decay mode and tagging category. The results of those checks are listed in Table [21 We observe 
no statistically significant asymmetry from fits to the control samples of non-CP decay modes. 

5 SUMMARY 

In summary, we report on improved measurements of sin2/3 and |A| that supersede our previous 
result §\. We measure sin2/3 = 0.710 ± 0.034(stat) ± 0.019(syst) and |A| = 0.932 ± 0.026(stat) ± 
0.017(syst). The updated value of sin2/3 is consistent with the current world average d6 and the 
theoretical estimates of the magnitudes of CKM matrix elements in the context of the SM |17j . 
The theoretical uncertainty on the interpretation of the measurement of sin2/3 in these modes is 
approximately 0.01 |18j . 
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